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Abstract 

The Cooper pair box qubit of a two-junction-SQUID configuration enables 
the readout of the qubit states by probing the effective Josephson inductance 
of the SQUID. This is realized by coupling the qubit to a high-Q tank circuit 
which induces a small alternating supercurrent in the SQUID loop. The effect 
of a small (but finite) geometrical inductance of the loop on the eigenstates of 
the system is figured out. The effect of qubit dephasing due to quadratic cou- 
pling to the tank circuit is evaluated. It is shown that the rate of dephasing 
in the vicinity of the magic points is relatively low unless the Josephson junc- 
tions forming the qubit are rather dissimilar. In the vicinity of the avoided 
level-crossing point such dephasing is always significant. 

1 Introduction 

The readout device is a critical component of any potential quantum computing 
circuit. For the Josephson qubits there is a number of sensitive cryogenic devices 
available (SQUIDs, switching Josephson junctions, single electron transistors and 
traps, etc.) enabling the readout of the qubit state. However, operation of these 
devices is usually associated with a significant exchange of energy between detector 
and qubit, so in order to avoid fast decoherence the detector must be reliably de- 
coupled from the qubit at the time of quantum manipulation. Recently, the class 
of Josephson qubit detectors based on the measurement of the reactive component 
of electrical signals related to nonlinear behavior of the Josephson inductance has 
been extensively studied Due to specific coupling to the qubit variables 
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Figure 1: Electric circuit diagram of the charge-phase qubit. The core of the qubit 
is the superconducting ring of small inductance L including two small tunnel junctions 
of capacitances Ci and C2 and Josephson coupling energies Eji and Ej2, respectively. 
The small island between the junctions is polarized by charge on small gate capacitance 
C g . The ring is inductively (mutual inductance M) coupled to an rf-driven parallel tank 
circuit. Voltage across the tank circuit is read out by a (presumably cold) preamplifier 
which is not shown. 



and non-dissipative characteristics of the Josephson supercurrent, these circuits 
can have a much weaker backaction and, therefore, cause lesser decoherence. More- 
over, these circuits may possibly enable quantum nondemolition measurements of 
a Josephson qubit [H]- 

In this paper, we consider the charge-phase qubit comprising a macroscopic 
superconducting ring including two small Josephson junctions with a small island 
in between (see Fig. 1), i.e., a Cooper pair box |S] of SQUID configuration. This 
setup is, in principle, similar to that of Quantronium 0, in which, however, a 
larger additional Josephson junction was inserted in the superconducting loop for 
the purpose of readout by switching this junction into a finite-voltage state. In 
our case, the ring of inductance L is inductively coupled to the tank circuit that 
enables both phase control and readout of the qubit. The special convenience of 
this simple setup is that it can function even without uncoupling the tank circuit 
from the qubit. It has recently been shown that operation of this qubit in the 
so-called magic points made it possible to decouple (in linear approximation) both 
the diagonal and the off-diagonal components of the variables of the measuring 
system (environment) from the qubit [H]. Moreover, the resonance shape of the 
noise spectrum allowed significant reducing relaxation and, in principle, performing 
single shot measurements. Here we analyze the quadratic effect of this noise and 
show that dephasing caused by this effect in the vicinity of the magic operation 
points is small, provided the qubit is sufficiently symmetric. 
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2 Qubit parameters and the model 

For the vanishingly small value of inductance L the behavior of the charge-phase 
qubit is primarily determined by the relationship between the energies associated 
with the island and the junctions. These are the Coulomb energy of the island, 
E c = e 2 /2Cs, and the Josephson coupling energies of the junctions, Ejij 2 = 
Wci,c2/2e, where Cs = C\ + C2 + C g is the island's total capacitance including 
the capacitances of the individual junctions C\ and C2 and the gate capacitance 
C g <C Ci,2; I c \,c2 are the nominal values of the Josephson critical currents. We 
assume that 

($ /27r)/ c0 = Ej = {E JX + E J2 )/2 ~ E c , (1) 

where $0 = h/2e is the flux quantum, Ejq and I c q = (I c ± + I C 2)/2 are the average 
Josephson coupling energy and the average critical current, respectively. The cou- 
pling energies of individual junctions are almost similar, Ej± ~ E.J2, so the values 
of the dimcnsionless parameters 312 = Eji j j2/2Ejo rs 0.5, while j± + 32 = 1. The 
junction capacitances are also assumed to be almost similar, C\ ~ C'2, and the 
dimensionless parameters = C^i/Cs ~ 0.5, while k\ + k 2 = 1. The qubit 
island is capacitively coupled to the gate voltage source V g which controls the po- 
larization charge Qq — C g V g . The inductively coupled coil Lt carrying current I 
induces the external magnetic flux $ e = ((f>o/27r)0 e = —MI applied to the qubit 
loop and produces the sum phase bias (f> = tpi + if2 — 4>e- The core Hamiltonian 
includes the charging (kinetic) energy and Josephson (periodic potential) term, 

H Q = AE c (n - Q /2e) 2 - Ej(<f>) cos x , (2) 

where the operator n of the number of excess Cooper pairs on the island is 

n = -i-^- = -i-^, [n,ip] = [n,x]=i- (3) 

The operators of the island phases are tp = (tpi — if 2)/ 2 and x — ¥' + 7( ( / ) )j where the 
angle deviation due to asymmetry of the junctions 7 = arctan [(ji — 32) tan(</>/2)]. 
The amplitude of the effective Josephson coupling energy of two junctions is equal 
to 

E.M) = (Ej, + Ej 2 + 2E.nE.j2 cos 0) 1/2 . (4) 

The two lowest eigenvalues of energy, Eq{Qq,4>) (ground state) and Ei(Qq,c/)) 
(first excited state), with corresponding eigenstates |0) and |1), form the basis 
suitable for qubit operation. In this basis the Hamiltonian Eq. (J2J takes the diagonal 
form 

Ho = -(e/2)<7„ (5) 

where energy e(Qo,4>) = foil = E± — Eq and a z is the Pauli matrix. The energy 
landscape is 2e periodic in Qq and 2n periodic in <j>. The plots can be found 
elsewhere (see, for example, Refs. jj] and |H]). The reverse Josephson inductance of 
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the qubit in the ground (excited) state L \ is determined by the local curvature of 
the corresponding energy surface, 

WQo>H(^J w . (6) 

Depending on the value of the ratio Ejq/E c and the relationship between the 
Josephson energies of the junctions, ji/j2, the inductance Lo.i can take either 
positive or negative values in different points on the Qo-4> plane. In the so-called 
magic points (extremum or saddle points), i.e., A (Qo = 0, 4> = tt)> B (Qo = 0, <j> = 
0) and C (Qo = e, = 0), the absolute values of Lq \ achieve local maxima, while 
in the avoided level-crossing point D (Qo = e, <f> = ir) its value is the largest [H]. 
Specifically, in point D, where effective coupling is small, Ej(tt) — \Ej\ — Ej%\ -C 
E c , these values for two states are equal to [5] 

The drive frequency u> r f (close to the resonance frequency of the tank circuit 
lot = (LtCt) 1 ^ 2 <C ksT/h) is much lower than the qubit frequency fi, so the 
induced classical oscillations of phase <fi are adiabatic. Due to coupling to the ring, 
the effective inductance seen by the tank circuit is ^U] 

L cS -±T- 7—7 L T~ J ■ (»J 

Li + ^0,1 -^0,1 

As a result, the resonance frequencies take the distinct values for the ground and 
excited states, w ,i = (L^^Ct)^ 1 ^ 2 , and this property is used for the radio fre- 
quency readout of the qubit pQ . 



3 Effect of finite inductance of the ring 

The problem of the inductance effects in the persistent current flux qubits had 
been addressed by Crankshaw and Orlando in Ref. where the corrections to 
the energy level values of these qubits were found. Recently, Maassen van den 
Brink had analyzed the finite inductance effect in the three-Josephson-junction 
jubits considering the self- flux as a "fast variable" ^21- He had shown that for 
small inductances this effect is merely reduced to renormalization of the individual 
Josephson coupling energies entering the system Hamiltonian. 

In the case of the charge-phase qubit, the sum phase (f> is also no longer a good 
variable. Instead, its combination with the normalized circulating current (which 
is a quantum-mechanical operator) is controlled by the external flux $ e , 

4>+^-LI s (^, X ) =4> e . (9) 
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The operator of the circulating current can be presented as [B] 

Is = I\\ 0) cos x + I± {4>) sin X- 



(10) 



In the given basis, the operator cosx is diagonal and sinx is off-diagonal, so the 
amplitudes of the longitudinal and transversal (Jj_) components of current I s 
are 



87T 
$0 



8tt . . . , E% 



(ji - h)(Kih + K ih) + 2jlj 2 (Kl - «2) cos 2 



(11) 



(12) 



Note that the operator I s is diagonal for arbitrary values of parameters Qo and 
only if the junctions are identical, i.e., j\ — 32 an d K\ = K2, so I± = 0. 
The additional magnetic energy term in the system Hamiltonian is 



H m = ($ /2tt) 2 



^ e ) 2 /2L = L/ 2 (0,x)/2. 



(13) 



The kinetic energy term associated with the charging of the chain of the qubit 
junctions connected in series, He = ( < &o/27r) 2 C , ^> 2 /2, where C — C qu bit + C stTay . 
Although the capacitance between the ends of the SQUID loop C s t r ay depends on 
the sample layout, it clearly dominates over the first term, C qu bit = C1C2/C5: which 
is typically about 1 fF The resonance frequency lul = (LC) 1 ^ 2 of the oscillator 
formed by L and C is therefore relatively low. Assuming Tiujl ~ ksT, but still 
u>l u) x f we recover classical behavior of <j>. In this case it is sufficient to take into 
account only the magnetic term associated with the circulating current given by 
Eq. P]l. 

The phase variable <f> is a single- valued function of 4> e for arbitrary values of Qo ; 
provided L is sufficiently small, viz., 

L\L^\<L\L^\(D)\ = -^-l3 L <l, (14) 
\3i - 32\ 

where the screening parameter is 

fa = ^LI c0 . (15) 

One can see that for an almost symmetric qubit (ji w j'2) the condition of smallness 
of inductance L, i.e. /3l < — J2I, is more stringent than that for ordinary rf 
SQUIDs, i.e. & <10. 

To derive the dependence (f>(4> e ), one can solve Eq. @ by iteration. The first 
step gives the expression 



0m&-(27t/*o)£I.(&,x)- 



(16) 
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Figure 2: Energies Eq and E\ calculated for the zero (wire frames) and nonzero (surface 
plots) values of screening parameter /3l for the case of mean Josephson coupling Ejo = 
+ Ej2) — 2E C and the asymmetry parameter a = ji — js = (Eji — Ej2)/{Ej 1 + 
E.J2) = 0.1. Black dots mark on the Qo-'&e plane the locations of the magic points A, B, 
and C and the avoided level-crossing point D. 



Then the corresponding one-dimensional Schrodinger equation, 

— - i-^-j - Ej(ct> e ) cosx- j [/||(0 e )cosx + /j_(</>e)sinx] 2 > ip = Eip, 

(17) 

is again of the type describing a quantum particle moving in a periodic potential. 
In the case of L = 0, this potential is harmonic and the equation is of the Mathieu 
type [TJ]. In the case of nonzero L, the potential energy includes a (small) second 
harmonic admixture. The eigenvalues and eigenfunctions of this equation are found 
by numerical methods. 

Figure 2 shows the two lowest eigenvalues of energy (the qubit levels) computed 
from the Schrodinger equation Eq. (|17|l . For the given parameters of the qubit and 
not very small value of (3l, the subsequent terms of iteration give substantially 
smaller correction. One can see that similar to the conclusion of Ref. JJ, the 
first order effect of inductance reduces to the lowering of both energy surfaces, 
Eq and E\. The correction to the level spacing Eq. 10) is, however, very small: 
e — > e + AeL, Ae^ <C e. As can be seen from Fig. 2, the effect of the deformation of 
the energy surfaces is most appreciable in the region of small effective Josephson 
coupling combined with charge degeneracy, e.g., in the vicinity of point D. Around 
the magic points A and - especially - B and C, which are most suitable for qubit 
operation, the change of local curvature of the surfaces is rather small. 
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4 Qubit dephasing 

The term in the system Hamiltonian which describes coupling of the qubit and the 
variable <5<I> e = MIf associated with the external flux noise induced in the ring by 
the tank circuit is given by 

#coupl = -Is (18) 

or, in the qubit basis Eq. JSJl, 

#coupl = -7^U z (J z -^Uy<7y = ~ [(en - C o)I|| {<t>e)cT z + 2s 01 Ij_(<f> e )<Jy] 5$ e , (19) 

where corresponding matrix elements are denoted by 

Coo = (0|cobx|0), c n = (l|cos X |l), and s 01 = |(0| sm X |l>|. (20) 

The remarkable property of the coupling described by Eq. Ijl9(l implies that the 
longitudinal term, U z oc In, vanishes in the optimal points A, B, C and D, i.e. 
where s'm(j> e = (see Eq. Ijllfl ') |16|. The transversal term, U y oc I±, in these points 
is, however, nonzero, and its low-frequency (<C £1) noise component can dephase the 
qubit due to renormalization of the level spacing ^2], e - * (e 2 + U 2 ) 1 ^ 2 ~ e + U 2 /2e, 
yielding the diagonal term 

- \xu z = -^a z = -^s 2 01 ll(ct> e )(6$ e ) 2 a z . (21) 

The problem of qubit dephasing due to quadratic longitudinal coupling to the en- 
vironment was addressed by Makhlin and Shnirman in Ref. II 7j where they focused 
on the experimentally relevant 1// and the Johnson- Nyquist noise power spectra. 
Here we extend this model to our case of a structured bath (low-frequency resonance 
tank circuit). 

For white noise // associated with losses in resistance Rt, the power spectrum 
density S$(ui) of variable 5<fr e is simply found from the network consideration. 
Neglecting small detuning of the tank due to coupling to the qubit ring (ujq ss oji ss 
lot) we arrive at the expression, 

S 9 (e>) = M 2 ^—— " T -oo < u < oo, (22) 

ttRt \oj z — U3 T Y + oj z oj t Q 1 

where Q = utCtRt = Rt/^tLt is the quality factor. Assuming that the noise is 
Gaussian we find the power spectrum density of fluctuations of the variable (SQ e ) 2 
(see, for example, Ref. [T%]). 

/+oo 
du>'S<s,(u')S<t,(w' -uj). (23) 
-oo 

Taking this integral at u> — > yields the value 



1 - r4^,3, . ( k BT 



2 



5 $2 (0) = -M 4 g^ T ( -jL- ) , (24) 
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which determines the power spectrum density Sx(0) of variable X, 

Sx(0) = ^ S 4 Q1 liS^(0), (25) 

and the rate of "pure" dephasing in the Bloch-Redfield approximation r v = 
Sx(0)/?i 2 . Finally, the qubit quality factor Q v = -k£1/T v can be presented as 



2 \soUj_J {E J0 k B Ty 

where k = M j '(LLt) 1 ^ 2 is the dimensionless coupling coefficient. 

To evaluate Q v for typical parameters of the qubit (leading for the operation 
points A, B and C to the value (c/Ejq) 3 ~ 10) let us assume that the quality factor 
of the tank circuit is Q ~ 300 while the product k 2 Q/3L ~ 10 which ensures sufficient 
resolution in determining the resonance frequencies luq.i. Taking the value of ratio 
E J0 huj T /(k B T) 2 ~ 0.1 (with E J0 = 100 /ieV, w t /2tt = 100 MHz and T = 1 K) we 
get an estimate: Q v ~ 10 2 (/ c o/soi^j_) 4 - 

In the optimal points A and D, the amplitude Jj_ (see Eq. I|12|l) is equal to 

(I±)a,d = 2(kij'i + K2ji)I c o ~ IcO, (27) 
while in points B and C 

{I±)b.c = 2 [(ji - J*2)(«iii + K2J2) + 2jij 2 (Ki - K2)] ^co 

~ (ii - J2 + Kl - K 2 )/ C 0- (28) 

The values of the matrix element sqi in points A and B are equal to [S] 

/ s _ \jiTh\Ej 

(soi)a,b « (29) 

while in points C and D 

(*oi)c,u « 0.5. (30) 

Summarizing, the order-of-magnitude estimates for parameter Q v in different 
operation points are 

(Q v ).4~10 4 (ji-j 2 )- 4 , (31) 

(Q<p)b ~ 10 4 (ji - 32 + ki - « 2 r 4 , (32) 

(Q ¥ ,)c - 10(ii - j 2 + ki - k 2 )~ 4 , (33) 

(Q v )i5 ~ 10. (34) 

For the Josephson junctions with a reasonably high symmetry in their parameters, 
say ±10% (e.g., j\ — j 2 = 0.1), the expected value of Q v in points A, B and C 
greatly exceeds the level 10 5 which should be sufficient for realization of active 
schemes for compensation of decoherence. In contrast to these figures, the value 
achieved in the avoided level-crossing point D is relatively low. 
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5 Conclusion 

We have shown that qubit dephasing due to nonlinear coupling in the vicinity of 
magic operation points is relatively weak if asymmetry of the junction parameters 
is sufficiently small and the effective temperature of the tank circuit sufficiently 
low. Since this temperature is associated with the backaction noise of a pream- 
plifier, its noise figure should be sufficiently low, as, for example, in the case of a 
cold semiconductor-based amplifier or a dc-SQUID-based amplifier [20]. Further 
improvement (if necessary) of the charge-phase qubit with radio frequency readout 
can be achieved by using a stage ensuring variable coupling between the ring and 
the tank circuit. This could be, for example, a flux transformer with a variable 
transfer function |21| . 

In conclusion, sufficiently symmetric parameters of the junctions make it pos- 
sible to effectively protect the charge-phase qubit against dephasing caused by its 
coupling to the radio- frequency readout circuit. The dominating mechanisms which 
may still cause appreciable dephasing in this kind of qubit are most likely related 
to intrinsic characteristics of the Josephson junctions (see, e.g., Refs. [22H2S])- 



6 Acknowledgment 

I would like to thank Alec Maassen van den Brink, Apostol Vourdas and Denis 
Vion for helpful comments. This work was partially supported by the EU through 
the SQUBIT-2 project. 



References 

[1] A. B. Zorin, "Cooper-pair qubit and Cooper-pair electrometer in one device," Physica 

C368, 284-288 (2002). 
[2] E. Il'ichev, A.Yu. Smirnov, M. Grajcar, A. Izmalkov, D. Born, N. Oukhanski, 
Th. Wagner, W. Krech, H.-G. Meyer, and A.M. Zagoskin, "Radio- frequency 
method for investigation of quantum properties of superconducting structures," 
cond-mat/0402559 and references therein. 
[3] A. Lupascu, C. J.M. Verwijs, R.N. Schouten, C.J. P.M. Harmans, and J.E. Mooij, 
"Nondestructive readout for a superconducting flux qubit," cond-mat/0311510 
[4] I. Siddiqi, R. Vijay, F. Pierre, CM. Wilson, M. Metcalfe, C. Rigetti, L. Frunzio, 
and M.H. Devoret, "An RF-driven Josephson bifurcation amplifier for quantum 
measurements," cond-mat/0312623 
[5] D. V. Averin, "Quantum nondemolition measurements of a qubit," Phys. Rev. Lett. 
88, 207901 (2002). 

[6] V. Bouchiat, D. Vion, P. Joyez, D. Esteve and M.H. Devoret, "Quantum coherence 
with a single Cooper pair," Phys. Scr. T76, 165-170 (1998). 

[7] D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. Urbina, D. Esteve and 
M.H. Devoret, "Manipulating the quantum state of an electrical current," Sci- 
ence 296, 886-889 (2002). 



10 



A. B. Zorin 



[8] A. B. Zorin, "Josephson charge- phase qubit with radio frequency readout: coupling 
and decoherence," Zh. Eksp. Teor. Fiz. 125(6), 1423-1435 (2004) [JETP 98(6), 
1250-1261 (2004)]. 

[9] W. Krech, M. Grajcar, D. Born, I. Zhilyaev, Th. Wagner, E. Il'ichev, and Ya. Green- 
berg, "Dynamic features of a charge qubit closed by a superconducting inductive 
ring," Phys. Lett. A 303, 352-357 (2002). 

[10] R. Rifkin and B. S. Deaver, Jr., "Charge-phase relation and phase-dependent con- 
ductance of superconducting point contacts from rf impedance measurements," 
Phys. Rev. B 13(9), 3894-3901 (1976). 

[11] D.S. Crankshaw and T. P. Orlando, "Inductance effects in the persistent current 
qubit," IEEE Appl. Supercond. 11(1), 1006-1009 (2001). 

[12] A. Maassen van den Brink, "Hamiltonian for coupled qubits," |cond-mat /0310425 

[13] In the rf-SQUID circuits comprising Josephson weak-links instead of tunnel junctions 
a sufficiently small value of effective capacitance C and, hence high ujl may cause 
quantum behavior of phase see, for example, R. Whiteman, V. Schollmann, 
M. Everitt, T.D. Clark, R.J. Prance, H. Prance, J. Diggins, G. Buckling, and 
J. F. Ralf, "Adiabatic modulation of a superconducting quantum interference 
device (SQUID) ring by an electromagnetic field," J. Phys.: Condens. Matter 
10, 9951-9968 (1998). 

[14] P. K. Hansma, Superconducting single-junction interferometers with small critical 

currents, J. Appl. Phys. 44(9), 4191-4194 (1973). 
[15] Handbook of Mathematical Functions, edited by M. Abramowitz and I. A. Stegun 

(U.S. GPO, Washingtin, D.C., 1972), Chapter 20. 
[16] In particular, quantum manipulation at a remarkably low dephasing rate of the 

Quantronium qubit was performed in operation point C 0. 
[17] Yu. Makhlin and A. Shnirman, "Dephasing of solid-state qubits at optimal points," 

cond-mat/0308297 

[18] S.A. Akhmanov, Yu. E. Dyakov and A. S. Chirkin, Vvedenie v Statisticheskuyu Ra- 
diofiziku i Optiku (Introduction to Statistical Radiophysics and Optics), (Nauka, 
Moscow, 1981), Chapter 5 - in Russian. 

[19] A. Shnirman, Yu. Makhlin and G. Schon, Noise and decoherence in quantum two- 
level systems, Phys. Scr. T102, 147-154 (2002). 

[20] M. Muck, J. B. Kycia and J. Clarke, "Superconducting quantum interference device 
as a near-quantum- limited amplifier at 0.5 GHz," Appl. Phys. Lett. 78(7), 967- 
969 (2001). 

[21] C. Cosmelli, M. G. Castellano, F. Chiarello, R. Leoni, D. Simeone, G. Torrioli, 
and P. Carelli, "Controllable flux coupling for the integration of flux qubits," 
cond-mat/0403690 

[22] E. Paladino, L. Faoro, G. Falci, and R. Fazio, "Decoherence and 1// noise in Joseph- 
son qubits," Phys. Rev. Lett. 88, 228304 (2002). 

[23] D. J. Van Harlingen, T. L. Robertson, B. L. T. Plourde, P. A. Reichardt, T. A. Crane, 
and J. Clarke, "Decoherence in Josephson-junction qubits due to critical current 
fluctuations," |cond- mat/0404307 



